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ABSTRACT 

We reconsider the emission properties of the BL Lac objects emitting in the high-energy 7- 
ray band exploiting the new information in the MeV-GeV band obtained by the Large Area 
Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope in its first three months 
of operation. To this aim we construct the spectral energy distribution of all the BL Lacs 
revealed by LAT and of the known TeV BL Lacs not detected by LAT, also including data 
from the Swift satellite, and model them with a simple one-zone leptonic model. The analysis 
shows that the BL Lacs detected by LAT (being or not already detected in the TeV band) 
share similar physical parameters. While some of the TeV BL Lacs not revealed by LAT have 
spectral energy distributions and physical parameters very similar to the LAT BL Lacs, a group 
of objects displays pecuUar properties (larger electron energies and smaller magnetic fields) 
suggesting different physical conditions in the emission region. Finally, we discuss possible 
criteria to effectively select good new candidates for the Cherenkov telescopes among the LAT 
sources, presenting a list of predicted fluxes in the very high-energy band calculated including 
the effect of the absorption by the extragalactic background light. 

Key words: radiation mechanisms: non-thermal — 7-rays: theory — 7-rays: observations - 
BL Lac objects: general 



' 1 INTRODUCTION 

The gamma-ray extragalactic sky at high (> 100 MeV) and very 
high (> 100 GeV) energies is dominated by blazars, associated to 
radio-loud active galactic nuclei with a relativistic jet closely ori- 
, ented to the Earth. The resulting relativistic amplification of the 
■ non-thermal jet emission makes blazars extreme objects, with ap- 
parent luminosities exceeding in the most powerful sources 10^*^ 
erg/s and variability timescales as short as minutes (e.g. Aharonian 
et al. 2007a, Albert et al. 2007a). 

Blazars are divided in two classes: i) flat spectrum radio 
quasars (FSRQs) show broad emission lines typical of quasars in 
their optical spectra and are the most powerful sources; ii) BL 
Lac objects, instead, display rather weak (or event absent) emis- 
sion lines (with, by definition, equivalent width less than SA) and 
have smaller luminosities. For both classes, the Spectral Energy 
Distribution (SED) is characterized by two broad peaks, generally 
interpreted as due to synchrotron (the low energy peak) and inverse 
Compton emission (at high energy; an hadronic origin for the high 
energy peak is supposed in the hadronic models, e.g. Muecke et al. 
2003). It is widely accepted that the SEDs of blazars follow a se- 
quence with the luminosity (Fossati et al. 1998, Donato et al. 2001): 
the powerful sources are characterized by SEDs whose synchrotron 
and IC peaks are in the submm/IR and MeV band, respectively. 



while the low-power BL Lacs show the maxima of their compo- 
nents in the optical/UV or even X-ray band and in the GeV- TeV 
region. Ghisellini et al. (1998) proposed that the sequence could be 
related to the balance between acceleration and cooling acting on 
the electrons emitting most of the power: electrons in the jet of FS- 
RQs are characterized by strong radiative losses, and thus cannot 
reach large energies. In BL Lacs, instead, the small cooling allows 
the electrons to be accelerated to very high energy, determiiring 
high frequencies for the synchrotron and IC emission components 
(see Ghisellini & Tavecchio 2008a, 2009 for a refined view). 

Blazars emitting in the TeV band are still a small group but 
their number is fastly increasing (see De Angelis et al. 2008 and 
Aharonian et al. 2008a for recent reviews^. The interest for these 
sources is driven by the possibility to get interesting clues on the 
acceleration processes of charged particles in relativistic flows and 
the possible use of their TeV emission to characterized the poorly 
known extragalactic background light (EBL; e.g. Stecker et al. 
1992, Stanev & Franceschini 1998, Mazin & Raue 2007). The fa- 
vorite candidates for the TeV detection are blazars peaking at the 
highest energies, called high frequency peaked BL Lacs (HBLs, 
Giommi & Padovani 1994), showing the peak of the synchrotron 
component in the UV-X-ray band and the IC peak close to 100 
GeV. For this reason most of the observation time has been de- 
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voted to observe HBLs (one of the most used list of candidates 
was presented in Costamante & Giiisellini 2002, hiereafter CG02). 
However, tfie new generation of Clierenkov telescopes now operat- 
ing (HESS, MAGIC, VERITAS) have enough sensibility to observe 
other, less bright in the TeV band, blazars. Indeed, a handful of IBL 
(intermediate BL Lacs) and LBLs (low frequency peak BL Lac ob- 
jects) belong to the present group ot TeV BL Lacs. One of the most 
critical points in the search for new TeV blazars is that the present 
Cherenkov telescopes have a quite limited field of view and cannot 
provide a full-sky survey. For this reason the selection of candidates 
has to rely on a variety of indirect methods. The simplest criterion 
(e.g. Stecker et al. 1996) is to choose the HBL whose synchrotron 
component peaks in the X-ray band assuming that the X-ray flux is 
a proxy for the flux of the TeV emission. A more refined method 
is based on two parameters, the X-ray and the radio flux (CG02). 
This method has been quite successful in helping to select the tar- 
gets, though not all the candidates have yet been detected. 

Recently, the Fermi collaboration released a catalogue of 
AGNs (with galactic latitude |6| > 10 deg) detected with high 
significance (lOcr above 100 MeV) by the Large Area Telescope 
(LAT) in the first three months of operation (Abdo et al. 2009a, A09 
hereafter). Excluding two radiogalaxies (Cen A and NGC 1275), 
the remaining 104 sources are all blazars (58 FSRQs, 42 BL Lacs 
and 4 with uncertain classificatior0). One interesting feature of this 
catalogue is the large number of BL Lacs. For comparison, the pre- 
vious EGRET list of BL Lacs detected in the 100 MeV - 10 GeV 
band (Hartman et al. 1999) comprised only 14 BL Lacs (out a to- 
tal of 60 blazars). The knowledge of the GeV part of the spectrum 
is rather important, allowing to better determine the shape and the 
position of the the peak of the emission of the TeV BL Lacs, when 
located between the bands covered by LAT and by the Cherenkov 
telescopes (see e.g. Aharonian et al. 2009). Moreover, the large 
number of BL Lacs, together with the good description of the spec- 
trum, is a useful tool to find new TeV candidates. 

The SEDs of the 23 most powerful blazars of the A09 list (with 
7-ray luminosity exceeding 10*** erg s^^) are presented and dis- 
cussed in Ghisellini, Tavecchio & Ghirlanda (2009). The emission 
models for the SEDs of all the A09 sources are presented in Ghis- 
ellini et al. (2009a). In this paper we focus our attention on the BL 
Lacs reported in the A09 list, including also the BL Lacs detected 
in the TeV band but not by LAT. We also discussed possible cri- 
teria to select, among the LAT BL Lacs, the best candidates to be 
detected by Cherenkov Telescopes. 

We use Ho = 70 km s"^ Mpc"\ Qa = 0.7, ST^m = 0.3. 



2 BL LAC OBJECTS IN THE FERMI/LAT 

The BL Lacs reported in A09 are reported in Table[T] where we also 
report whether there are Swift observations, whether the source was 
detected by EGRET (or there is an upper limit, Fichtel et al. 1994) 
and whether the source has been detected in the TeV band. We ex- 
clude from this list 6 sources (AO 0235+164, PKS 0332-403, PKS 
0426-380, PKS 0537-441, PKS 1057-79, S5 1803-^784), classified 
as BL Lac in A09, but showing broad emission lines with luminosi- 
ties typical of quasars (see Ghisellini et al. 2009a for details). More- 
over, for three sources reported in Table[T](0141+268, 1054+2210 



^ Note that one of the sources classified as FSRQ is a radio-loud Narrow 
Line Seyfert 1 galaxy (Abdo et al. 2009b, Foschini et al. 2009). 
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Table 1. The Fermi BL Lac objects in the A09 list. We have excluded 6 
sources, defined as BL Lacs based on the small equivalent width of the 
lines, but showing broad emission lines typical of quasars (see Ghisellini et 
al. 2009a). In the last 3 columns we indicate if there are Swift observations; 
if the source was detected by EGRET (UL stands for an upper limit given by 
Fichtel et al. 1994) and if it is a TeV source, u: redshift uncertain; S: lower 
limit from Sbai'ufatti et al. (2005); n: not studied in this paper due to lack 
of multiwavelength data; a: redshift from NED. A09 reports z = 0.416. 

and 1250+532) the paucity of the data prevents even a rough de- 
scription of the SED and we excluded them in the following analy- 
sis. Therefore the total number of LAT BL Lacs considered here is 
33. This list contains 12 BL Lacs detected in the TeV band, exactly 
half of the known TeV BL Lacfl. 

A09 report the spectral parameters (flux and slope) obtained 
with a power law model fit to the average LAT spectrum obtained 
during the first three months. The spectral characteristics of the TeV 
(black triangles) and non-TeV BL Lacs (red circles) are compared 
in Fig.[T] in which we report the photon index and the flux above 
100 MeV resulting from the power law fit. Clearly, TeV sources 
preferentially lay on the lower part of the plot, indicating a harder 

^ Note, however, that lES 2344+514 lies at a galactic latitude |6| < 10 
deg, not covered by the A09 list. 
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Table 2. The known TeV emitting BL Lacs not in the A09 list. * |fe| < 10 
deg. 



GeV spectmm than that of the non TeV ones. This is not surpris- 
ing, since most of the BL Lacs detected so far at TeV energies are 
HBLs with the high-energy peals; of the SED above 100 GeV, im- 
plying a hard spectmm in the LAT band. Only BL Lac itself and S5 
0716+71 (and, marginally, 3C66A and W Comae) have a soft pho- 
ton index. Occasionally, BL Lac displays a very hard GeV spectrum 
(e.g. Nandikotkur et al. 2007). The hardest (and faintest) source is 
MS 1050.7+4946 {z = 0.140), withaphoton index T^, = 1.4±0.2. 
TeV and non-TeV sources share the same distribution of fluxes. 
This fact ensures that the difference in the distributions in slope is 
an intrinsic feature and not the result of a selection effect related 
to the higher sensitivity of LAT for harder sources (dotted line in 
Fig[T] see also A09). 

In Table |2] we report the 12 TeV BL Lacs not present in the 
A09 list. Neither of these sources was detected by EGRET. In the 
following we will also consider these sources for the modeling of 
the SED. 



3 ANALYSIS OF THE SWIFT DATA 

Several of the BL Lac in A09 (and of the non-LAT TeV sources) 
have been observed by Swift. 

Swift/XKT (Burrows et al. 2005) and Swift/VWOT (Roming et 
al. 2005) data have been reduced using the standard procedure (de- 
scribed in, e.g., Ghisellini et al. 2007), with the latest version of 
the HEASOFT (6.6.3) and CALDB (2009/5/6) packagefl XRT data 
have been fitted with XSPEC12.In most cases a single power law 
model absorbed with a column density fixed to the Galactic value 
(taken from Kalberla et al. 2005) provides a good fit of the data. 
In few cases a broken power law (again with a Galactic column 
density) is required. Results are shown in Table g] UVOT magni- 
tudes have been corrected for Galactic absorption using the values 
of Schlegel et al. (1998) for the V, B and U filters and the formulae 
by Pei (1992) for the UV filters and converted into fluxes following 
Poole et al. (2008). Results are reported in Table|5] 



4 SPECTRAL ENERGY DISTRIBUTIONS 

We assembled the SEDs of the BL Lacs in the A09 and of the 12 
TeV BL Lac not included there using historical (mainly from the 



http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ 



2.4 



2.2 



1.8 



1.6 



1.4 



D03a , 

. ' 060 ; 




tl216 I 
"••^ T ^155 

T421 



J I I I I I 



_l I I I I I I I 



_l I I I J 



1 10 
FiooMev [10-« Ph cm-^ s-'] 



Figure 1. Gamma-ray photon index versus the integral flux above 100 MeV 
for the BL Lac objects detected by Fermi/LAI (A09). Black triangles refer 
to BL Lacs detected in the TeV band (see Table [T). TeV detected sources 
belong to the hardest sources, while the distribution of fluxes is very similar 
for both groups. The dashed line reports the limit sensitivity for the sources 
reported in A09. Note that the LAT sensitivity increases with the hardness 
of the spectra. 



NASA Extragalactic Database^ and the BZCAT/ASDS archivtjfl), 
FermifLKI from A09 (when available) and, if available. Swift data 
taken in (or close to) the period of the Fermi observations (August 4 
- October 30, 2008). When a large amount of data are available for 
a given source (for instance the well studied Mkn 421, Mkn 501 and 
PKS 2155-304) we choose to plot the most representative datasets, 
to show the possible range of variations in the different bands. For 
clarity, we report only the TeV spectra already corrected for the 
effect of the absorption through the interaction with the EBL. Ab- 
sorption is taken into account considering only one of the mod- 
els discussed in literature, the "lowSFR" model of Kneiske et al. 
(2004), based on a low estimate of the EBL close to values recently 
derived (e.g. Aharonian et al. 2006). 

For the sources detected by LAT we plot the result of the fit 
with a power law model to the average spectrum (lower bow tie) us- 
ing the parameters reported in A09. To show the variability we also 
report the spectrum (upper bow-tie) obtained by using the weekly 
averaged peak flux provided by A09 and assuming the same slope 
of the average spectrum. This last assumption is probably violated 
in some of the sources, since the spectrum likely changes with the 
flux (possibly in a rather complex fashion, see e.g. Aharonian et al. 
2009). As already discussed, the LAT spectra for HBLs are gen- 
erally hard, falling onto the raising part of the high-energy peak. 
In IBL/LBL sources the LAT spectra rather probe the peak of the 
emission. We warn the reader that, as also discussed in A09, not 
always the power law model provides a good fit of the data over 
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the entire LAT band. In particular, the slope is often dictated by the 
low energy bins, characterized by the smallest errors. In general, we 
expect that the GeV spectra of HBLs soften approaching the peak 
of the bump. A clear example is given in Aharonian et al. (2009) 
for PKS 2155-304. For IBLs and LBLs the situation could be even 
more complex, since the LAT band is likely centered on the high- 
energy peak. For the TeV BL Lacs not present in A09 we report an 
upper limit in the GeV band calculated using the sensitivity curve 
reported in A09 (dotted line in FiglTJ. 

For the sources detected by LAT it is interesting to compare 
the LAT and the TeV spectra. In general, for the sources with a 
good coverage at TeV energies (especially Mkn 421, Mkn 501, PKS 
2155-304, lES I959-I-650), the LAT spectrum agrees very well with 
TeV spectra taken at low level. For the sources with a limited num- 
ber of observations/detections in the TeV band, the LAT data seems 
to connect rather smoothly with the available TeV spectrum. 

Finally, another interesting point to consider concerns the vari- 
ability in the GeV band. The great majority of these BL Lac are 
known to be extremely variable in the TeV band, showing varia- 
tions of factors as large as 10 on timescales as short as hours or 
even minutes (e.g. Aharonian et al. 2007, Albert et al. 2007a). On 
the other hand, variability in the LAT band seems to be more mod- 
erate (less than a factor of 2 in flux), as shown by the peak flux 
spectra. However, variability on timescale shorter than one week is 
probably diluted and thus not visible in these data. 

Our aim is mainly to reproduce "average" SEDs of the 
sources, to investigate these BL Lacs as a population. For this rea- 
son, in the case of the sources detected by LAT we choose to calcu- 
late the model using the LAT spectrum as a guideline to reproduce 
the high-energy peak. Indeed, the 3-months average LAT spectrum 
is a more reliable indicator of the average state than the fastly- 
evolving TeV spectrum. This is also reinforced by the fact that 
published TeV spectra are often biased, being taken during high 
levels of the emission. LAT data, instead, are more representative 
of the average state of the sources. As noted above, in most cases 
the model is not very distant from the observed TeV spectra at low 
levels. In these cases we try to locate the position of the high-energy 
peak between the LAT and the (soft) TeV spectrum. For the sources 
not detected by LAT we model the SED using the TeV spectrum. 

When possible we fix the synchrotron component by using 
Swift data (providing a good coverage of the crucial optical-UV 
and X-ray band) taken during the three months covered by the LAT 
observations (or close to the TeV observations for the sources not 
detected by LAT). When Swift data exactly covering the period 
of high-energy observations are not available we choose to use, if 
available, the data closest in time. In most cases they have been ob- 
tained few months after or before the observations of LAT or TeV 
telescopes. However, we are aware that, given the extreme vari- 
ability displayed by these sources, especially at X-ray frequencies, 
these SEDs should be taken with some cautions. For instance, for 
two cases (PKS 2155-304, 3C66A) there are more than one Swift 
observations in the period covered by the LAT showing important 
variations of the X-ray emission, both in flux and slope. In these 
cases we calculated two models, corresponding to the extremes of 
the X-ray flux. For PG 1553-F113 and RGB J0152-I-017 there are 
no Swift observations close in time with LAT. However, we find 
archival observations showing huge variations in the X-rays. For 
PG 1553-1-1 13, for which no X-ray observation nearly simultaneous 
with the 7-ray data exits, we produced two different models using 
the two extreme X-ray spectra. For RGB JO 152-1-0 17 we use the X- 
ray spectrum (taken with XMM and RXTE, Aharonian et al. 2008) 
close in time with the TeV observations. In the case of MrklSO, 



for which we have three Swift observations (one in spring 2008, 
two just after the LAT observing period) showing huge variability 
we model the low level X-ray spectrum and the MAGIC spectrum. 
Finally, for the interesting case of 1426-1-428 we present two mod- 
els, one referring to the average HEGRA and CAT spectrum (Aha- 
ronian et al. 2003a) and one reproducing the recent simultaneous 
Swift and MAGIC data (Leonardo et al. 2009). 



5 MODELING OF THE SEDS 

We use the one-zone synchrotron self-Compton model fully de- 
scribed in Maraschi & Tavecchio (2003). Briefly, the emitting re- 
gion is a sphere with radius R with a tangled and uniform magnetic 
field B. The relativistic effects are described by the Doppler factor 
S. The (purely phenomenological) distribution of the emitting rel- 
ativistic electrons is described by a broken power law model with 
normalization K and indices ni from 7miii to 7b and 712 above 
the break up to 7max. The model includes the full Klein-Nishina 
cross-section for the calculation of the inverse Compton spectrum, 
particularly important for the emission above the GeV band (e.g. 
Tavecchio et al. 1998). 

We recall that with one-zone models, such as that adopted 
here, one can not reproduce the emission at the longest wave- 
lengths, since the emission is self absorbed below the millime- 
ter band. That part of the SED is due to outer regions of the jet, 
not important for the modeling of the high-energy emission. Multi 
zone emission models, such as those discussed in the framework 
of structured jet models (Georganopoulos & Kazanas 2003, Ghis- 
ellini et al. 2005), could be required for the detailed modeling of 
the SED even at 7-ray energies (e.g. Anderhub et al. 2009). How- 
ever, with the available non-simultaneous data the one-zone model 
can satisfactorily reproduce most of the SEDs. 

Models are shown in Figs. I6l7l and[8l and the parameters used 
to reproduce the data are reported in Table |6] We assume a red- 
shift of 2 = 0.3 for all the sources with unknown distance. For the 
sources with more than one model we indicated the parameters for 
the low (0 and high (h) state. 

The parameters of this simple one-zone model can be com- 
pletely constrained when the value of peak frequency and flux of 
both emission components and an estimate of the minimum vari- 
ability timescale are available (Tavecchio et al. 1998). In most cases 
the synchrotron peak can be constrained by the optical and X-ray 
data, while the IC peak is less determined, since the LAT spectra 
generally cover only the raising (or in some cases the declining) 
part of the bump. In the fit procedure we decide to locate the IC 
peak just above (or below) the LAT band. With this choice (see 
also below) we minimize the required Doppler factor. 

Inspection of Table|6]shows that the parameters are quite sim- 
ilar in most of the sources, with some notable exceptions. In Fig.|2] 
we summarize the results showing the values of the magnetic field 
and the electron break Lorentz factor 7b. Open circles (blue trian- 
gles) report the values corresponding to TeV BL Lacs detected (not 
detected) by LAT, while green open squares indicates LAT BL Lacs 
not detected in the TeV band. Most of the LAT sources are local- 
ized in the region corresponding to magnetic field between 0.1 and 
1 G and 7b is between 10'' and 10^. BL Lac itself requires a slightly 
larger magnetic field and a lower value of the break Lorentz factor. 
The reason is that its peaks (both synchrotron and IC) are located at 
lower energies (note that it is characterized by one of the steepest 
LAT spectrum, Fig[T| locating the IC peak well below 100 MeV). 

Sources not detected by LAT seem to clearly separate, being 
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Figure 2. Magnetic field versus the break Lorentz factor 7b for the BL Lacs 
modeled with the one-zone synchi'otron-SSC model. Red open circles (blue 
filled triangles) show the values for the known TeV sources detected (not 
detected) by LAT. The two populations are clearly divided, with the sources 
not detected by LAT populating the region with small magnetic field and 
large electron Lorentz factor. Open green squares are for the LAT BL Lacs 
not detected in the TeV band. 



characterized by low magnetic field and large electron Lorentz fac- 
tor. Inspection of the SEDs reveals that these sources could be di- 
vided in two groups: (i) sources rather similar to those detected 
by LAT, whose parameters fall quite close to the region popu- 
lated by the LAT BL Lacs; (ii) peculiar sources: in this group 
there are particular sources (lES 0229-1-200, lES 1101-232, RGB 
J0152-I-017, lES 0347-121, H1426-I-428 and, even if not so ex- 
treme, H 2356-309), characterized by a very hard TeV continuum 
(de-absorbed). The interpretation of these spectra is still a debated 
issue (e.g. Katarzynski et al. 2005, Sleeker et al. 2007, Aharonian 
et al. 2008b, Boettcher et al. 2008). To reproduce the SED of these 
sources one needs to use rather extreme parameters, in particular a 
rather large energy for the emitting electrons. A possible model in- 
vokes the emission from relativistic electrons distributed as a power 
law in energy with a very high value of the minimum energy. In 
these conditions it is possible to get a rather hard IC continuum, 
F{v) oc v^^^ (see Katarzynski et al. 2005 and Tavecchio et al. 
2009 for a detailed discussion of the cases of lES 1101-232 and 
lES 0229-1-200, respectively). 

Based on these differences, we can predict that sources be- 
longing to the group (i) will probably be detected in deeper obser- 
vations of LAT. Instead, sources of type (ii) are intrinsically differ- 
ent: the presence of the very hard IC continuum, peaking above 1 
TeV, makes very difficult for LAT to detect them, even with pro- 
longed exposures. In this sense, LAT is probably able to probe 
only a portion of the TeV BL Lac population, i.e. the sources with 
the peak located not far from 100 GeV. All sources with a peak at 
higher energies will probably escape the detection of LAT. This fact 
has important consequences for the use of LAT data in the search 
for new TeV candidates. 

In Fig.[3](upper panel) we show the distribution of the Doppler 
factors used in the models (the shaded histogram is for non-TeV 
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Figure 3. Distributions of the Doppler factors (upper panel) and minimum 
vaiiability timescales (lower panel) derived through the modelling of the 
SED. The shaded area indicates the distribution of sources not detected in 
the TeV band. 



sources). Most sources have 5 between 20 and 30. Few extreme 
sources require 5 > 40. TeV sources requiring these extreme values 
are mainly those with a very hard TeV spectrum. Large 5 in the case 
of non-TeV sources are required for those BL Lacs with hard LAT 
spectra and synchrotron peak in the optical band (see below). 

In the bottom panel of Fig.[3]we show the minimum variability 
timescale (as measured in the observer frame), fvar = (1 + z)R/c5 
derived from the modeling. The bulk of the sources have tvar 
around 1 and 10 hours. These are typical timescales predicted with 
the standard one-zone SSC model. Faster variations (such as those 
observed in PKS 2155-304 and Mkn 501, Aharonian et al. 2007, 
Albert et al. 2007a) require modifications of this scheme (e.g. Ghis- 
ellini & Tavecchio 2008b, Ghisellini et al. 2009b, Giannios et al. 
2009, Neronov et al. 2008). 



5.1 Comments on specific sources 

Inspection of SEDs reveals few cases (e.g. 1717-1-177, 2322-1-396) in 
which the model does not provide a good fit of the LAT spectrum. 
All these sources are characterized by a synchrotron peak located 
close to the optical band and a hard LAT spectrum, implying a IC 
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peak above 10 GeV. In these conditions, the SSC model requires 
a rather large Doppler factor, 5 > 50, and a low magnetic field 
(e.g. Tavecchio & Ghisellini 2008 for a discussion). In the case of 
2322+396 even with 5 = 40 (and B = 0.007 G; corresponding 
to the lowest open square in Fig. |2] ) the fit is not completely sat- 
isfactory. A possibility to solve this problem is to assume that, as 
already discussed, the bow tie is somewhat misleading, since the 
power law is not a good model for the high-energy bins of the ac- 
tual LAT spectrum. Therefore, contrary to the impression given by 
the power law fit, the peak would be located at few GeV, alleviating 
the requirements on the model parameters. 

Notable is the case of 1514-241 (AP Lib). In this case it is 
rather difficult to reconcile the extremely steep optical-UV con- 
tinuum, the hard X-ray spectrum and the (non simultaneous) LAT 
spectrum. The model strongly underpredicts the X-ray flux. The 
first obvious solution to this problem is to assume that the non si- 
multaneous X-ray spectrum is not representative of the real state 
of the X-ray emission during the LAT observations. Another possi- 
bility to reconcile the X-ray and the LAT spectra is to assume that 
the emission originates in a structured jet (Ghisellini et al. 2005), in 
which the radiative interplay between a fast inner core (the spine) 
and a slower outer layer would result in a complex spectrum, with 
more than one component contributing in the X-ray-7-ray band. 
The detailed discussion of this point is beyond the scope of this 
paper. 

Another interesting source is S5 0716-F71. Although the SSC 
model reproduces well the SED when only the LAT data are consid- 
ered, the fit of the TeV emission observed in April 2008 by MAGIC 
although feasible (Anderhub et al. 2009), is difficult, requiring a 
rather broad high-energy component, extending above 300 GeV. 
Also in this case a viable solution invokes the emission from a 
structured jet (Tavecchio & Ghisellini 2009, Anderhub et al. 2009). 



6 PROSPECTS FOR THE DETECTION IN THE TEV 
BAND 

The search for new TeV BL Lacs is a rather difficult task for 
Cherenkov telescopes. Their limited field of view does not allow an 
all-sky survey. The method usually applied is then to select good 
candidates based on indirect estimators of the level of the TeV 
emission, such as the X-ray luminosity (e.g. Stecker et al. 1996). 
With the launch of Fermi and the availability of an all-sky survey at 
7-rays up to hundreds of GeV we have now the unprecedented pos- 
sibility to choose the best TeV candidates selecting those sources 
showing an intense emission up to 10-50 GeV (the maximum en- 
ergy probed by the LAT with moderate exposures). This selection 
based on the LAT flux and spectrum is the most effective way to 
select TeV candidates, since it uses information on the level and 
spectrum of the emission just below the threshold of the current 
Cherenkov telescopes. Fig. [T] could thus already provide some in- 
teresting clues on some of the best TeV candidates. Apart BL Lac 
and S5 0716-1-71 (characterized by relatively soft LAT spectra), the 
other known TeV BL Lacs appear to be bright and relatively hard in 
the LAT band. Within this group, 1215-^303, 1717-^177, 2322-1-396 
and Ap Lib, should be considered as optimal candidates. An out- 
standing source is MS 1050.7-1-4946. Its exceptionally hard spec- 
trum makes it an optimal candidate but the flux is rather low, prob- 
ably limiting the possibility to detect it. 

Fig. |4] shows the LAT photon index as a function of the X- 
ray flux. Here the TeV sources (again, excluding BL Lac and 
S50716-I-71) are located in the region of hard GeV spectra and 
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Figure 4. Photon index of the LAT spectrum versus the X-ray flux at 1 keV 
(in vFi,) for the LAT BL Lacs. Black triangles indicate the BL Lac detected 
in the TeV band). Besides the two IBL BL Lac and S507 16-1-71 afl the TeV 
BL Lacs, characterized by a hard spectrum (F < 2) and a bright X-ray 
continuum, are confined in the area defined by the dashed fine. 



large X-ray fluxes. This is simply due to the fact that the ma- 
jority of the known TeV sources (mostly HBLs) are bright X-ray 
sources and have the IC peak above 100 GeV, thus showing a hard 
GeV spectrum. Notably, some objects not yet detected in the TeV 
band lye in the same locus of the TeV BL Lacs, that we can ap- 
proximately define by the rectangle plotted with dashed lines. To- 
gether with some of the sources already noted above (I2I5-I-303, 
I7I7-I-I77, 2322-1-396), there are other sources that should be con- 
sidered as ideal TeV candidates: 0447^39, 0502-1-675, 0133-1-388, 
0033-1938. Note also that 1424-1-240, very recently detected by 
VERITAS and MAGIC (Ong et al. 2009, Teshima et al. 2009), 
belongs to this group of sources. This reinforce the idea that the 
sources included in the area defined by the dashed lines should be 
considered the best candidates for the TeV detection. Note again 
the peculiar position of MS I050.7-I-4946. 

Finally, it is also interesting to compare the positions of the 
LAT BL Lacs with the other blazars in the radio-X-ray flux plot 
of CG02. They proposed to select the TeV candidates among the 
blazars with the largest X-ray and radio fluxes. The theoretical idea 
behind this choice is that the level of the synchrotron X-ray emis- 
sion is a proxy for the number of high-energy electrons, while the 
radio flux is a measure of the energy density of the soft photons 
providing the seeds for the IC process. The brightest TeV BL Lacs 
should be those with the largest radio and X-ray fluxes. This meth- 
ods have been proved to be quite effective: indeed most (though not 
all) the TeV BL Lacs detected in the past years belong to the ob- 
jects selected by CG02. In Fig.[5]we report a revision of the CG02 
plot including all the LAT BL Lacs for which both radio and X-ray 
fluxes are available. TeV BL Lacs (shown by triangles) are concen- 
trated in the left-top side of the diagram. Quite interestingly, almost 
all the LAT BL Lacs are characterized by a large radio flux. We in- 
terpret this fact as a confirmation of the idea behind the CG02 plot: 
the brightest 7-ray emitters are in fact those with the largest radio 
flux. The large span in the X-ray flux, instead, comes from the pres- 
ence of both HBL and LBL in the A09 list, conventionally sepa- 
rated in this plane by the line at aiix = 0.75 (Giommi & Padovani 
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Figure 5. X-ray versus radio flux for the BL Lacs in the A09 list. Filled 
triangles indicates known TeV BL Lacs. Gray points are from the original 
plot of Costamante & Ghisellini (2002). Dashed lines indicate the X-ray 
flux level of 1/iJy (horizontal) and the radio flux of 1 Jy (vertical). The 
diagonal dashed line indicates the locus characterized by an effective slope 
between the radio and the X-ray band of ajjx = 0.75, conventionally 
separating HBLs (ajjx < 0.75) andLBLs (ctjix > 0.75). 



1994). Also in this diagram we can extract some particularly inter- 
esting candidates looking at the area occupied by the known TeV 
BL Lacs: again we find 1215-1-303, 0447-439, 0133-1-388. Two other 
outstanding sources are 0033-0192 and 0502-1-675. They are the 
sources with the lowest radio flux among the LAT BL Lacs, but 
they have rather hard GeV spectra and large X-ray flux (see Fig|4ll, 
suggesting an important emission at TeV energies. 

An important caveat in considering the results of these plots is 
that, when selecting the best candidates, one must also consider the 
effect of the absorption by the EBL, especially important above 50 
GeV and for sources at redshift z > 0.1. It is possible that sources 
appearing good candidates accordingly to the plots discussed above 
are actually penalized by strong absorption due to a relatively large 
distance. 

In order to provide a criterion to select sources for TeV tele- 
scopes including also the effect of absorption, we calculated the ex- 
pected fluxes of all the sources in the band covered by Cherenkov 
telescopes by using two independent way to estimate the intrinsic 
flux of the source: (i) in the first case we simply assume that the 
spectrum measured by LAT extends unbroken up to TeV energies; 
(ii) as a second estimate we use the high-energy spectrum predicted 
through the application of the emission model to the SEDs (in Fig. 
[8] we show both the intrinsic and the absorbed spectrum). Of course 
the first choice is rather generous, since we know that the actual 
spectrum of TeV blazars generally steepens above 100 GeV. On the 
other hand, in the second case the predicted fluxes are rather small, 
since the model IC component generally decreases quite fast above 
the LAT band due to the decreasing scattering efficiency. However, 
as discussed above, we cannot exclude the contribution from the 
emission of other components at TeV energies, thus resulting in 
larger fluxes than those expected with our simple one-zone mod- 
elling. Considering both estimates should thus give a reasonable 
idea of the real TeV fluxes from these sources. Also here absorption 
is calculated using the "lowSFR" model of Kneiske et al. (2004). 
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Table 3. Predicted flux above 50 and 300 GeV extrapolating the LAT spec- 
trum (L) and from the SED model (A/) for the LAT BL Lac not yet detected 
in the TeV band. Fluxes above 50 GeV are in units of 10~^^ ph cm"^ s~^, 
fluxes above 300 GeV in units of 10^^'' ph cm""^ s^^. For sources with 
unknown redshift we use z = 0.3. 

The results are reported in Table[3] We report the fluxes above 
50 GeV and 300 GeV predicted with the extrapolation of the LAT 
spectrum (L) and from the model (Af). We recall that for all the 
sources with unknown redshift we fix it to 2: — 0.3. With the largest 
threshold only few sources could be accessible with current instru- 
ments, while lowering the energy at 50 GeV all the sources should 
be characterized by fluxes larger than 10^^^ ph cm^"^ s^^. Not 
surprisingly, the favorite sources are those already noted above, 
in particular 0447-439, 0502+675, 1215+303 and 1717+177. MS 
1050.7+4946 is the preferred source accordingly to both ways to 
estimate the intrinsic flux. Less bright, but still interesting are 0033- 
1938, 0133+388, 2322+396. In particular the first source, though 
showing a rather intense intrinsic TeV emission, is located at large 
redshift (z = 0.61), determining a relatively low observed flux. 



7 SUMMARY AND CONCLUSIONS 

We have studied all the BL Lac objects detected by the LAT instru- 
ments onboard Fermi in its first three months of operations. Among 
them there are 12 sources detected also in the TeV band. To these 
we also added in our study the remaining 12 TeV BL Lacs not de- 
tected by LAT. For all these sources we have built the SED using 
also optical-UV and X-ray data from Swift. The SEDs have been 
modeled with a simple one-zone leptonic model in order to derive 
the main physical parameters of the emission region. We found that 
the sources have rather similar Doppler factors, distributed around 
5 = 2Q ~ 30. Some exceptional source require larger values, up to 
5 — 40. Typical magnetic fields are in the range B = 0.1 — 1 G, 
while the Lorentz factor of the electrons emitting most of the power 
(at the peaks of the SED) are between 7b = 10^ and 10^. 

Interestingly, while a fraction of TeV BL Lacs not yet detected 
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by LAT share the same values of the parameters with the LAT de- 
tected objects, several sources are characterized by rather extreme 
values (very small B and large electron Lorentz factor), suggesting 
perhaps a real difference in the physical conditions of the emission 
region in these sources. In fact, different possibilities to produce 
hard spectra have been discussed (see Tavecchio et al. 2009 for a 
discussion). In particular, Katarzyfiski et al. (2005) showed that if 
the electrons (assumed to follow a power law distribution with en- 
ergy) are characterized by a relatively large value of the minimum 
energy, a very hard spectrum can be achieved (F^ cx v^^^). The 
same hard spectrum would be visible in the synchrotron compo- 
nent below the soft X-rays. Stecker et al. (2007), instead, pointed 
out that, under specific conditions, diffusive shock acceleration pro- 
duces power law electron distribution much harder than the canon- 
ical 7~^. In both cases there would be the possibility to have SSC 
spectra with slopes harder then the canonical value V — 1.5, the 
limiting slope commonly assumed in deriving upper limits to the 
EBL (e.g. Aharonian et al. 2006). 

The full-sky survey at 7-rays performed by LAT offers a 
unique opportunity for the study of the high-energy emission of 
BL Lacs. For sources already discovered in the TeV band, we can 
now study in detail the high-energy component in the GeV band, 
and how it connects with the emission at higher energies. The com- 
parison between the variability in the two bands (GeV and TeV) is 
expected to give important clues on the emission mechanisms and 
dynamics of the relativistic electrons in the jet (e.g. Aharonian et 
al. 2009). The possibility to have a continuous monitoring of the 
sources can provide good alerts for the Cherenkov telescopes, al- 
lowing a better understanding of the global behaviour of the high 
energy emission of BL Lacs including their duty cycles. 

Particularly interesting, as we have shown in this paper, is the 
prospect to use the LAT data to select good candidates to be ob- 
served by Cherenkov telescopes. The first three months of observa- 
tions already provided several good sources worth to be considered. 
Due to the strong variability of these objects at high-energies, is 
likely that many other interesting BL Lacs, even with bright emis- 
sion, would be revealed in the future. However, the LAT selection is 
probably biased against the sources with extremely hard TeV spec- 
tra, for which our model predict a rather low flux in the MeV-GeV 
band, well below the capabilities of LAT. 
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02/03/2009 


2.3 ±0.4 


- 


- 


8.6/4 


0.08 


1101-H384 


3/12/2008 


1.7 ±0.5 


0.8 ±0.22 


2.3 ±0.2 


23/140 


96.2 


iziy-t-zoj 


//D/zUUS 


z.4o it U.UO 










1424-1-240 


11-19/6/2009 


2.2 ± 0.08 


1.24 ±0.21 


3.0 ±0.1 


413/491 


3.1 


1514-241 


08/04/2007 


1.7 ± 0.3 






3/12 


0.65 


15429-H6129 


18-20/1/2009 


2.5 ± 0.4 






0.3/2 


0.05 


1553-Hll 


6/10/2005 


2.24 ± 0.03 






203/221 


72.8 




20/4/2005 


2.38 ± 0.06 






84/73 


27.0 


1652-1-398 


12/5/2008 


1.68 ± 0.15 


0.9 ± 0.14 


2.1 ± 0.06 


129/126 


16.4 


1717-H177 


08/01/2009 


1.72 ± 0.37 






_* 


0.16 


1749-H096 


26/02/2007 


2.13 ± 0.21 






23/50 


17.6 


1959-H650 


27/10/2008 


2.07 ± 0.05 






46/48 


86.0 


2005-489 


8/10/2007 


2.69 ± 0.14 






13/11 


5.7 


2136-428 


19/12/2008 


2.3 ±0.4 


- 


- 


2.6/2 


0.05 


2155—304 


17/8/2008 


2.63 ± 0.04 






100/91 


85.9 




17/10/2008 


2.39 ± 0.08 






42/33 


37.9 


2200-H420 


29/8/2008 


1.97 ±0.1 






34.6/41 


6.4 


2322-H396 


2//3/2009 








* 


0.05 


0152-1-017 


2/12/2007 


2.23 ± 0.12 






10/14 


5.2 




13/1/2008 


2.33 ± 0.18 






8/11 


0.48 


0229-H200 


8/8/2008 


1.8 ±0.1 






23.6/21 


10.3 


0347-121 


3/10/2006 


1.95 ±0.05 






80.4/70 


4.4 


0548-322 


28/11/2006 


1.7 ±0.1 


1.8 ±0.7 


2.1 ±0.3 


78/71 


21.0 


0710-1-591 


26/2/2009 


1.80 ± 0.05 






70.52 


47.0 


0806-H524 


13/3/2008 


2.5 ±0.1 






16/16 


9.26 


1133-H704 


10/5/2008 


1.7± 0.2 


1.2 ±0.2 


2.4 ±0.2 


40/36 


17.8 




30/10/2008 


2.44 ±0.1 






12.5/14 


8.4 




24/11/2008 


1.35 ±0.5 


1.1 ±0.3 


2.0 ±0.1 


91/92 


52.3 


1426-H428 


9/6/2008 


2.51 ±0.23 






212/219 


5.1 



Table 4. Results of the X-ray analysis. [1]: name. [2]: observation date (dd/mm/yyyy). [3]: photon index. [4]: break energy (broken power law model). [5]: 
high energy photon index (broken power law model). [5]: value of the and degrees of freedom. [6]: unabsorbed 0.3-10 keV flux in units of 10~^^ erg 
cm~2 s~^. * poorly determined spectrum, the C-Statistics was used. 
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source 


OBS date 
dd/mm/yyyy 


V 


B 


U 


UVWl 


UVM2 


UVW2 


00311-1938 


21/1 1/2008 


16.38 ± 0.03 


16.70 ± 0.02 


15.80 ±0.02 


15.70 ± 0.02 


15.67 ±0.02 


15.74 ±0.01 


0048-097 


04/06/2008 


— 


— 


— 


— 


— 


15.16 ±0.01 


0109-H22 


31/05/2006 


14.79 ± 0.01 


15.23 ± 0.01 


14.38 ± 0.01 


14.46 ± 0.01 


14.45 ± 0.01 


14.56 ±0.01 


0219-1-428 


2/10/2008 


14.46 ± 0.07 


14.73 ± 0.01 


14.73 ± 0.01 


13.97 ± 0.01 


— 


14.06 ±0.01 




3/10/2008 


14.39 ± 0.01 


14.74 ± 0.01 


13.89 ± 0.01 


13.99 ± 0.01 


14.00 ± 0.01 


14.11 ± 0.01 


0301-243* 


24/3/2009 


15.46 ± 0.03 


15.78 ± 0.02 


14.91 ± 0.02 


14.88 ± 0.03 


14.77 ± 0.04 


14.89 ± 0.03 


0447-439 


19/12/2008 


14.36 ± 0.01 


14.62 ± 0.01 


13.69 ± 0.01 


13.56 ± 0.01 


13.43 ± 0.01 


13.51 ± 0.01 


0502-H675 


04/01/2009 


16.6 ±0.01 


17.01 ± 0.01 


16.20 ±0.01 


16.35 ± 0.01 


16.46 ±0.01 


16.49 ± 0.01 


0712-H5033 


21//1/2009 0.218 


17.04 ± 0.04 


17.47 ± 0.03 


16.77 ± 0.03 


16.96 ± 0.04 


17.09 ± 0.05 


17.16 ±0.04 


0716-1-714 


28/11/2008 


— 


— 


— 


— 


— 


17.71 ± 0.03 


0735-H178 


18/05/2007 


— 


— 


— 


— 


16.00±0.01 


— 


0814-H425 


04/01/2008 


— 


— 


— 


— 


— 


18.27 ±0.02 


0851-H202 


11/11/2008 


14.85 ± 0.02 


15.32 ±0.02 


14.61 ±0.02 


14.75 ± 0.02 


14.79 ± 0.02 


15.01 ±0.02 


1011+496 


8/5/2008 


15.28 ±0.03 


15.53 ±0.02 


14.58 ± 0.01 


14.42 ± 0.01 


14.23 ± 0.02 


14.30 ± 0.01 


1050.7+4946 


02/03/2009 


— 


— 


21.74 


17.53 ± 0.05 


— 


— 


1101+384 


3/12/2008 




— 


— 








1219+285 


7/6/2008 


14.68 ±0.01 


15.05 ± 0.01 


14.20 ±0.01 


14.17 ±0.01 


14.01 ± 0.01 


14.15 ±0.01 


1424+240 


11-19/6/2009 


14.49 ± 0.02 


14.81 ±0.01 


14.49 ± 0.02 


13.97 ±0.01 


13.9 ±0.01 


14.03 ± 0.01 


1514-241 


08/04/2007 


— 


16.02 ± 0.03 


15.67 ±0.02 


16.15 ±0.01 


— 


16.39 ± 0.09 


15429+6129 


18-20/1/2009 


16.40 ± 0.04 


16.73 ± 0.03 


15.87 ±0.03 


15.88 ± 0.04 


15.83 ± 0.04 


15.95 ±0.03 


1553+11 


6/10/2005 


13.94 ± 0.01 


14.28 ±0.01 


13.34 ± 0.01 


13.34 ± 0.01 


13.29 ± 0.01 


13.42 ± 0.01 




20/4/2005 


14.26 ±0.01 


14.68 ±0.01 


13.69 ±0.01 


13.74 ±0.01 


13.72 ± 0.01 


14.01 ±0.01 


1652+398 


12/5/2008 


13.94 ± 0.01 


14.66 ± 0.01 


14.16 ± 0.01 


14.16 ± 0.01 


14.06 ± 0.01 


14.12 ± 0.01 


1717+177 


08/01/2009 


17.71 ± 0.09 


18.01 ± 0.05 


17.27 ± 0.04 


17.36 ± 0.04 


17.37 ± 0.05 


17.59 ± 0.03 


1749+096 


26/02/2007 


14.19 ± 0.01 


14.76 ± 0.01 


14.11 ± 0.01 


14.60 ± 0.01 


14.98 ± 0.02 


15.06 ± 0.01 


1959+650 


27/10/2008 


15.09 ± 0.03 


15.54 ±0.02 


14.84 ±0.02 


15.04 ± 0.03 


15.12 ± 0.03 


15.14 ±0.02 


2005-489 


8/10/2007 


13.32 ± 0.01 


13.71 ± 0.01 


12.75 ± 0.01 


12.66 ± 0.01 


12.54 ±0.01 


12.67 ±0.01 


2136-428 


19/12/2008 


15.42 ± 0.03 


15.77 ±0.02 


14.90 ±0.02 


14.92 ± 0.03 


14.86 ± 0.04 


14.98 ± 0.03 


2155-304 


17/8/2008 


12.87 ±0.01 


13.18 ±0.01 


12.20 ±0.01 


12.07 ±0.01 


11.91 ±0.01 


12.01 ±0.01 




17/10/2008 


13.44 ± 0.01 


13.75 ± 0.01 


12.79 ±0.01 


12.69 ±0.01 


12.56 ±0.01 


12.67 ±0.01 


2200+420 


29/8/2008 


14.87 ±0.01 


15.69 ±0.01 


15.37 ±0.01 


16.04 ± 0.02 


16.72 ± 0.03 


16.73 ± 0.02 


2322+396 


2/3/2009 


18.7 ±0.3 


18.9 ±0.1 


18.5 ±0.1 


18.8 ±0.1 


19.1 ±0.1 


19.6 ±0.1 


0152+017 


13/1/2008 


15.59 ±0.03 


16.38 ± 0.03 


15.84 ± 0.03 


15.94 ± 0.03 


15.86 ± 0.03 


15.94 ± 0.02 


0229+200 


8/8/2008 


17.04 ±0.07 


18.11 ±0.07 


17.92 ± 0.08 


18.21 ± 0.08 


18.45 ±0.1 


18.38 ± 0.06 


0347-121 


3/10/2006 


17.18 ± 0.07 


17.65 ± 0.05 


16.86 ± 0.04 


16.54 ±0.04 






0548-322 


28/1 1/2006 


16.3 ± 0.03 


16.96 ± 0.02 


16.61 ± 0.03 


16.72 ± 0.03 


16.68 ± 0.04 




0710+591 


26/02/2009 


16.56 ± 0.05 


17.21 ± 0.04 


16.48 ± 0.04 


16.41 ± 0.04 


16.34 ± 0.05 


16.40 ±0.03 


0806+524 


13/03/2008 


15.78 ± 0.04 


16.17 ±0.02 


15.29 ± 0.02 


15.27 ±0.02 


15.29 ± 0.02 


15.32 ±0.02 


1133+704 


30/10/2008 










15.41 ± 0.03 


15.42 ±0.02 


1426+428 


09/06/2008 


16.39 ± 0.04 


16.99 ± 0.03 


16.21 ±0.01 


15.99 ±0.02 


15.89 ±0.03 


15.86 ±0.02 



Table 5. Results of the UVOT analysis. (*) Average of two observations performed the same day. 
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source 


7min 


7b 


O'niax 


ni 


"2 


B 


K 


R 


S 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


UZiy+42f5 / 




1 /I \y 1 a4 
1.4 X iU 


Q \y 1 A^J 
O X iU 


z 




i.l 


Q \y 1 a4 
o X iO 





2j 


/I 




o o 1 a4 
2.2 X iU 


Q \y 1 a5 
o X lU 


o 
z 


/I Q 


U.o 


\y 1 a4 

X io 


D 


ZJ 


0716+71 


1 


6 X 10^ 


8 X 10^ 


2 


4.5 


0.2 


o o s/ 1 a4 
2.2 X lU 


8.3 


47 


1 Al 1 1/1 OA 


30 


8 X lO'^ 


2 X 10^ 


1.9 


4.3 


0.33 


-1 , , -1 r\3 
l.D X lU 


8.5 


25 


1 1 A 1 1 1 /I 


1 


10^ 


4 X 10"^ 


2 


4 


0.24 


O.O X lU 


4 


20 


1 o 1 n 1 oo ^ 
iziy+zoj 




1.5 X IC^ 


1 a6 
10 


2 


4.2 


0.3 


4 X 10'^ 


3 


25 


1424+240 


3 X 10^ 


2 X 10'^ 


2 X 10^ 


2 


AA 


0.55 


1.3 X 10'^ 


7.6 


35 


1553+113 / 


1 


4 X lO'^ 


3.5 X 10^ 


2 


4.2 


0.4 


2 X 10^ 


7.5 


20 


/i 


1 


3 X lO'^ 


5 X 10^ 


2 


3.8 


2 


3.35 X 10^ 


9 


23 


1652+398 


1 


1.2 X 10^ 


6 X 10^ 


2 


3.9 


0.313 


4.5 X 10^ 


1 


20 


1959+650 


1 


5.7 X 10-^ 


6 X 10^ 


1.9 


3.4 


0.4 


7 X 10^ 


7.3 


18 


2005-489 


10^ 


1.3 X 10-^ 


10^ 


2 


4.8 


0.7 


1.6 X 10^ 


8 


22 


2155-304 / 


5 X IQi 


1.5 X lO'^ 


10*^ 


2 


4.5 


0.28 


1.65 X 10-* 


17.5 


15 


/i 


5 X 10^ 


2.5 X lO'^ 


10^ 


2 


4.5 


0.9 


1.5 X 10^ 


14.5 


20 


2200+420 


10 


9 X 10^ 


2 X 10^ 


2 


3.9 


1.4 


8 X 10^ 


2 


15 


Ul j2+01 / 


10 


3 X lO'' 


iU 


2 


3.5 


0.03 


10 


4 


25 


r\nA 1 "tr\n 

0229+200 


7.5 X 10^ 


5 X 10*^ 


4 X 10^ 


2.3 


3.3 


5 X 10^'^ 


y.5 X 10 


35 


40 


AQ /I "7 1 T 1 

UJ4 /— 121 


3 X 10'^ 


2 X 10^ 


4 X 10^ 


2.2 


3.5 


0.0025 


2.2 X 10 


60 


44 


r\c /I o Tio * 

05 4 o— 322 


1.5 X 10^ 


2.8 X 10'"^ 


1.5 X 10^ 


2 


4.1 


0.1 


10 


3.6 


25 


A"? 1 A 1 ^A 1 * 

U / lO+jyl 


2.5 X 10^ 


5 X 10"* 


10 


2 


4. 1 


0.13 


2.5 X 10*^ 


4 


25 


0806+524 


4 X 10^ 


5 X lO"* 


2 X 10^ 


2 


4.5 


0.39 


8 X 10^ 


5.7 


20 


1101-232 


4 X 10"* 


3.5 X lO""' 


2 X 10^ 


1.8 


4.3 


0.02 


1.9 X 10^ 


13 


41 


1133+704 


lO'' 


5 X lO"* 


2 X 10^ 


2 


4 


0.15 


1.2 X 10* 


4 


20 


1218+304 


10* 


1.3 X 10^ 


10^ 


2 


4.5 


0.2 


1.2 X 10^5 


1 


35 


1426+428 / 


102 


6 X lO* 


10*5 


2 


4 


0.18 


2.8 X 10^ 


5.9 


25 


h 


7 X 10-^ 


1 X lO"* 


6 X 10*^ 


2 


3.1 


0.033 


5 X 10^ 


3 


35 


2344+514 


1 


10" 


7 X 10^ 


2 


3.2 


0.1 


3 X 10* 


3.5 


25 


2356-309 


10^ 


2 X 10^ 


3 X 10'^ 


2 


3.8 


0.12 


6 X 10^ 


3 


35 



Table 6. Input model parameters for all the TeV BL Lacs with (upper panel) and without (lower panel) LAT detection. [1]: source. [2], [3] and [4]: minimum, 
break and maximum electron Lorentz factor. [5] and [6]: slope of the electron energy distribution below and above 7;,. [7]: magnetic field [G]. [8]: normaHzation 
of the electron distribution in units of cm"''. [9]: radius of the emission zone in units of 10^^ cm. [10]: Doppler factor. *: no high energy data yet available. 



source 


Tmin 


7b 




ni 


"2 


B 


K 


R 


5 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


0033-1938 


100 


2 X 10"' 


10^* 


2 


4 


0.4 


2 X lO"' 


6.5 


26 


0048-097 


100 


5 X 10^ 


2 X 10^ 


2 


4.2 


0.46 


9 X 10"' 


7 


18 


0109+22 


30 


7 X 10^ 


7 X 10"' 


2 


4.5 


0.2 


3.5 X 10"' 


7.6 


30 


0118-272 


10 




1.5 X 10= 


2 


4.6 


0.35 


2 X 10"' 


10 


25 


0133+388 


1 


2 X lO* 


10^ 


2 


4.2 


0.8 


5.5 X 10^ 


8 


20 


0301-243 


10 


lo-* 


2.5 X 10^ 


2 


4.6 


0.33 


1.4 X 10"' 


10 


2 


0447^39 


100 


1.5 X lO"* 


10'5 


2 


4.4 


0.4 


2 X 10"' 


5.1 


20 


0502+675 


5 X 10^ 


7 X lO"* 


10'^ 


2 


4.1 


0.75 


3.5 X 10^ 


10 


20 


0712+5033 


60 


lo-* 


3 X 10"' 


2 


4.5 


0.17 


1.3 X 10^ 




25 


0735+178 


30 


lo-* 


8 X 10"' 


2 


4.6 


0.13 


7 X 10"' 


8 


25 


0814+425 


70 


2 X lO* 


1.5 X 10^ 


2 


4.6 


0.042 


3.3 X 105 


6 


25 


0851+202 


250 


2 X 10^ 


5 X 10"' 


2 


4.2 


0.5 


4 X 10^ 


7 


18 


1050.7+4946 


7 X 10^ 


10^ 


5 X 10** 


2 


4.6 


0.025 


10^ 


3.8 


20 


10586+5628 


30 


8 X 10^ 


10« 


2 


4 


0.1 


7 X 10"' 


5 


25 


1215+303 


10^ 


2 X 10" 


10^ 


2 


4.4 


0.1 


9 X 10"* 


4.5 


20 


1514-241 


1 


2 X lO"* 


5 X 10"' 


2 


4.9 


0.012 


10* 


10 


40 


15429+6129 


40 


1.4 X lO"* 


3 X 10^ 


2 


4.8 


0.085 


2.3 X 10* 


9 


28 


1717+177 


500 


8 X 10^ 


6 X 10^ 


2 


3.6 


0.05 


6.3 X 10^ 


3 


21 


1749+096 


200 


2 X 10^ 


1.05 X 10^ 


2 


4 


0.75 


1.6 X 10^ 


6 


25 


2136-428 


200 


1.4 X 10"' 


9 X 10* 


2 


4.8 


0.12 


3 X 10* 


9 


28 


2322+396 


50 


3 X 10"' 


9 X 10* 


2 


4.9 


0.007 


8 X 10* 


9 


40 



Table 7. Input model parameters for the LAT BL Lacs not yet detected in the TeV band. See Table|6]for definitions. 
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Figure 6. Spectral Energy Distribution of the BL Lacs detected by LAT and TeV telescopes. Historical data are from NED. Optical-UV and X-ray data from 
Swift are from this work. Bow ties report the power law model fitting the LAT data (Abdo et al. 2009). See text for more details. TeV spectra have been 
deabsorbed with the "lowSFR" model of Kneiske et al. (2004). References for the TeV data: 0219+428: Acciari et al. (2009a); 0716+71: Anderhub et al. 
(2009); 1011+496: Albert et al. (2007b); 1101+384: Zweerink et al. (1997), Maraschi et al. (1999), Albert et al. (2007c); 1219+285: Acciari et al. (2008); 
1424+240: the datapoint is the observed (not deabsorbed) flux reported by Ong et al. (2009); 1553+113; Albert et al. (2007d); 1652+398; Albert et al. (2007a). 
The solid line is the result of the one-zone leptonic model. 
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Figure 6. -continue- 1959+650: Tagliaferri et al. (2008), Aharonian et al. (2003b); 2005-489: Aharonian et al. (2005a); 2155-304: Aharonian et al. (2005b), 
Aharonian et al. (2007a); BL Lac: Albert et al. (2007e). 
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Figure 7. Spectral Energy Distribution of tlie BL Lacs detected by TeV telescopes but not by LAT. TeV spectra have been deabsorbed with the "lowSFR" 
model of Kneiske et al. (2004). References for the TeV data: 0152+017: Aharonian et al. (2008c); 0229+200: Aharonian et al. (2007b); 0347-121: Aharonian 
et al. (2007c); 0548-322 and 0710+591: no TeV spectrum pubhshed yet; 0806+524: Acciari et al. (2009b); 1101-232: Aharonian et al. (2006); 1133+704: 
Albert et al. (2006a). The solid line is the result of the one-zone leptonic model. 
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Figure 7. -continue- 1218+304: Albert et al. (2006b); 1426+428: Aharonian et al. (2003), Leonardo et al. (2009); 2344+514: Albert et al. (2007f); 2356-309: 
Aharonian et al. (2006). 
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Figure 8. Spectral Energy Distribution of tlie BL Lacs detected by LAT but not by TeV telescopes. The solid line is the result of the one-zone leptonic model, 
the dashed line is the same model taking into account the absorption of 7-rays by the extragalactic background light (using the Kneiske et al. 2004 lowSFR 
model). 
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Figure 8. -continue- 
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